Electromagnetic form factors of the nucleon in the chiral quark soliton model 
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In this paper we present the derivation as well as the numerical results for the electromagnetic 
form factors of the nucleon within the chiral quark soliton model in the semiclassical quantization 
scheme. The model is based on semibosonized SU(2) Nambu - Jona-Lasinio lagrangean, where the 
boson fields are treated as classical ones. Other observables, namely the nucleon mean squared 
radii, the magnetic moments, and the nucleon-A splitting are calculated as well. The calculations 
have been done taking into account the quark sea polarization effects. The final results, including 
rotational 1/N C corrections, are compared with the existing experimental data, and they are found 
to be in a good agreement for the constituent quark mass of about 420 MeV. The only exception is 
, the neutron electric form factor which is overestimated. 
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I. INTRODUCTION 



In the last years one of the most challenging problems in elementary particle physics seems to be the solution of 
QCD in the low energy region. The main difficulties are due to the non-perturbative effects caused by the growing 
effective coupling constant of the fundamental theory in the low energy limit. This prevents one from using the 
well-known main tool of theoretical physics — the perturbation theory. Because of this, the most intriguing features 
of QCD, confinement and chiral symmetry breaking, still remain conceptual and practical problems. 

The above mentioned obstacles have initiated an increasing interest in non-perturbative methods and effective 
low-energy models of hadrons. The effective models are expected to mimic the behavior of QCD at low energies. 

The simplest purely fermionic Lorentz invariant model describing spontaneous chiral symmetry breakdown is the 
Nambu-Jona-Lasinio (NJL) model 0. The NJL lagrangean contains chirally invariant local four-fermion interaction. 
In its simplest SU(2) form it has the following structure: 

Cnjl =V(ip-m )y + ^. + (W 75 *)2 ] t (1) 

where ^ is the quark field, G is the coupling constant, f are the Pauli matrices in the isospin space and too is the 
current mass the same for both up and down quarks. 

Applying the well known bosonization procedure Q the NJL model is expressed in terms of auxiliary meson fields 
a, 7?: 

a 2 

* [ ip- mo - 9(a + i 7? • ] * - + 7T 2 ). (2) 

Here, g is the physical pion-quark coupling constant implying that 7? are the physical pion fields. We assume that 
the meson fields are constrained to the chiral circle: 

o\x)+t\x) = fl , (3) 

where = 93 MeV is the pion decay constant. 

In fact, an equivalent effective chiral quark meson theory || can be derived from the instanton model of the QCD 
vacuum. 

In the chiral quark soliton model @]5||| based on the lagrangean (||) (frequently referred to simply as NJL model) the 
baryons appear as a bound state of N c valence quarks, coupled to the polarized Dirac sea. Actually, the baryon solution 
of the model is obtained in two steps. In the first step, motivated by the large N c limit, a static localized solution 
(soliton) is found by solving the corresponding equations of motion in an iterative self-consistent procedure |^|3j, 
assuming that the a and tt fields have hedgehog structure. However, this hedgehog soliton does not preserve the 
spin and isospin. In order to describe the nucleon properties one needs nucleon states with proper spin and isospin 
numbers. To this end, making use of the rotational zero modes, the soliton is quantized ^,0. Within the semiclassical 
quantization scheme quite successful calculations have been reported for the nucleon-delta mass splitting || and axial- 
vector form factor |j| as well as some results for the nucleon electric form factors |l0|. V ery recently, in the semiclassical 
quantization procedure important 1 /N c rotational corrections have been derived |llf| for isovector magnetic moment 
as well as for the axial coupling constant. These corrections are derived in a formalism (TTJ , based on the path 
integral approach Q. The scheme involves a time-ordering of collective operators which follows from the collective 
path integral. It improves considerably the theoretical values for the isovector magnetic moment, which is strongly 
underestimated ]T^ | in the leading order. 

Besides its numerical success, there are also theoretical reasons why the NJL model, described by lagrangean (|^), 
is considered as the one of the most promising effective theories describing low energy QCD phenomena. 

First, the model is the simplest quark model describing spontaneous breaking of the chiral symmetry, a basic feature 
of QCD. This allows to treat the pions as Goldstone particles and to include familiar methods of the current algebra 
into the model. As for the confinement, which is absent in the model, some recent lattice QCD results for the vacuum 
correlation functions of hadronic currents |Dj ] suggest that the confinement mechanism seems to have little impact 
on the light hadron structure, in particular on the nucleon and Delta structure. 

Second, there are various hints P, p||l5| ] that the NJL-type lagrangean can be obtained from QCD in various low- 
energy approximations. It should be stressed that the large N c limit plays a prominent role in those considerations. 

The aim of this paper is to calculate the electromagnetic nucleon form factors within the SU(2) chiral quark soliton 
model, based on the semibosonized NJL-type lagrangean, where the vacuum polarization effects taken into account. 
The calculations done in the semiclassical quantization scheme will include the rotational l/N c corrections, which 
have been shown pT| to be important for the isovector magnetic moment. In the derivation of the matrix element 
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of the electromagnetic current in the semiclassical quantization scheme not-commuting collective operators appear 
and the final result depends on the ordering of these operators. In fact, similar to ref. Jll]] the problem of the 
proper ordering of the collective operators is solved in the present paper using a path-integral approach M to the 
semiclassical quantization, which allows a natural access to the problem. This scheme allows to calculate the nucleon 
electromagnetic form factors, as well as such quantities like the electric mean squared radii, the magnetic moments, 
the nuclcon-A energy splitting, and the electric and magnetic charge distributions of the nucleon. 

The paper is organized as follows. In Section 2 we derive general expressions for the nucleon form factors in 
the NJL model. They are split in Dirac sea and valence quark contributions both including rotational corrections, 
coming from the expansion of the effective action in powers of the angular velocity of the soliton. Since the NJL 
model is non-renormalizable an ultraviolet regularization, which preserves the electromagnetic gauge invariance, is 
introduced in Section 3. The general formalism developed in Sections 2 and 3 is then applied to evaluate the nucleon 
electromagnetic form factors in Section 4. In Section 5 we present and discuss our numerical results. 



II. DIRAC SEA AND VALENCE QUARK CONTRIBUTIONS TO NUCLEON FORM FACTORS 

Our main goal in this section is to derive a general formula for nucleon form factors in the NJL model including 
rotational corrections. We start with the path integral representation for the nucleon matrix element of the corre- 
sponding quark current. After integrating out the quark fields we find that the result splits in two parts, the Dirac 
sea and the valence quark contributions. Usually the remaining integral over the meson fields is evaluated within 
the saddle point approximation. However, the saddle-point classical solution has a hedgehog symmetry that breaks 
both rotational and SU(2) isospin symmetries. In order to restore them one has to consider a rotating soliton and to 
perform a functional integration over the time-dependent orientation matrix of the soliton. Expanding in powers of 
the angular velocity up to linear order we derive rotational corrections to the nucleon form factors. A special attention 
is paid to the correct ordering of non-commuting collective operators describing the rotational degrees of freedom of 
the soliton. In this section the consideration is limited to the non-regularized case and the ultraviolet regularization 
is postponed till Section 3. 

Under the chiral circle restriction (|^) we can introduce the SU (2) chiral meson field 

U = -^((7 + in ■ t) (4) 

and work with the Euclidean lagrangean 

L = ¥D(U)^, (5) 

where the operator 

D(U) = 74(7^ + ML/ 75 + m ) = d T + h(U) (6) 

includes the Euclidean time derivative d T (r = X4) and the Dirac one-particle hamiltonian reads: 

h(U) = 74 (7fcd fc + MU^ + m ). (7) 

Here M = gf n is the constituent quark mass. For convenience, we choose to work with hermitean Euclidean Dirac 
matrices 7^ = j a . We use notation 

v n = l±Jl v + IzT*tft. (8) 
2 2 w 

The aim of this paper is to calculate the electromagnetic nucleon form factors. However, to be more general we start 
with a matrix element of quark current ^O^l , where O is some matrix with spin and isospin indices, and represent 
this matrix element by the Euclidean functional integral [Q with the lagrangean (^|): 

(N'J\¥(0)CW(0)\N,p)= lim \ [ d 3 xd 3 ye~ tifs ' +l P s 

t^+oo A J 

Jw Jv^l J P* t J A r,(x',-T/2)* t (0)O*(0) J^(f,T/2)e"/ d4z * tD(c/) *. (9) 
Here the nucleon state is created by the nucleon current J N (x) constructed of N c (number of colors) quark fields \& 
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Jn{x) = ^ e c -- c ^r { ^; T ^ } vi/ clQ1 (x)...vi/ cjVcajVc (x) (10) 

where Cj are color indices and r"j 3 1S a matrix with cti standing for both flavor and spin indices. J and T denote 
the nucleon spin and isospin, respectively. In the limit of large Euclidean times x' A ,Xi only the lowest nucleon state 
survives in (|). 

We assume that the constant Z in eq. (^) is chosen so that the nucleon states obey the non-relativistic normalization 
condition: 

(N,j/\N,p) = (2t:) 3 S^(j/ -p). (11) 
Let us integrate the quarks out in (^) . The result is naturally split in valence and Dirac sea parts (see Fig.|l|) : 

(N 1 \i?\^OV\N,p) = (N',p , \^0^\N,p) sea + (N',p'\^0^\N,p) val , (12) 



where 



-ip x +ipx 



(N\0\^O^\N : p) sea = N CT \im i Jd 3 xd 3 ye 
JvUT%--^T%^*\{(x',Tl2\^ ) \x,-Tl2) 0kak Sp{o 7 y (0, 0| — 10, 0)y 7 } 



fe=i 

xe N c Trlo S [D(U)] (13) 



and 



(N'^l^O^iO^NtP)™ 1 = N c lim i / d 3 a;d 3 j/e-^ /+ipt 

t^+oo Z J 

r Nc 1 



x(x\T/2|^|0,0) ft7 O 7y (0,0|^|f,-T/2) yQ1 e^ Trl °^)]. (14) 

In the limit of large number of colors N c the functional integrals over the chiral field U(x) in cqs. (|l3|), ( |l4|) can be 
evaluated by the saddle point method. The saddle point solution U is a static localized field configuration with the 
hedgehog symmetry 

U[x) = e iP(\*\X* a T*)/\A , (15) 

where P(|x|) is the profile function of the soliton. This saddle point solution breaks the space rotational and isospin 
symmetries. Therefore, even in the leading order in N c , one should go beyond the field U, extending the path integral 
over all fields of the form: 

U{x) = R{x A )U{x)R\x A ) , (16) 

where R(xi) is the time-dependent SU(2) orientation matrix of the soliton. 
For this ansatz (p^ ) the operator D(U) (g) can be written as 

D(U) = R[D(U) + R t R]R j < . (17) 

The dot stands for the derivative with respect to the Euclidean time £4. Similarly, the quark propagator in the 
background meson field U can be rewritten as 

(x I — !— I x') = R(xA (x I r I x') RUxi) . (18) 

\ 1 D{u) v 4M 1 D(U) + R^R ' ' V 47 K ' 

Let us insert ansatz ([!(]) into our formulas for the matrix element (|l3|), (jlj) and replace the functional integral over 
U(x) by the integral over the time dependent orientation matrix R(xi): 
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(N',p'\^0-$>\N,tf) sea = N Crr \im i J d 3 xd 3 ye- lifs ' +l P s 



VRT%~P N °r%- aN °* n [i?(T/2)(f ,T/2\ 



k=l 



D(U) 



\x-T/2)tf{-T/2) 



Sp{o 7 y ^(0X0,01—1-10,0)^(0)] ( } 



£>(t/)' 



, e N c Ti'\og[D(U)+RiR] 



7 7 



(19) 



and 



(A", ^| (0) | JV, ^> 



val 



N c lim ^ / d^d 3 ^-^'^" 



• / Pfle^TrioglD^+fltfl]^...^^...^* ^Qj i? ( T/ 2)(f' jr /2|_l^|f,-r/2)i? t (-r/2) 



fc=2 



R(T/2)(x', T/2\^-\0, 0)flt(O)] O^y [(0, Ol^yla?, -T/2)R){-T/2) 



l'otl 



(20) 



In eqs. (|19|), (J20J) despite of ansatz ( |16[ ) we still have a complicated path integral over VR. In order to solve this 
problem we make essentially use that in the large N c limit the angular velocity R'R of the soliton as well as its 
derivative is suppressed Q. Hence, similarly to ref. Q we can use expansions in powers of it 



N c Tr log[D(U) +R t R] = N c Trlog[L>(£/)] + <d sea J drSp(R^R) 2 



and 



1 rtR- 1 



D{U) + R!R D(U) D(U) D(U) 



(21) 



(22) 



In eq. (|T]) the second term of the action is a function of the angular velocity only. As we will see later, the coefficient 
Q sea is the Dirac sea contribution to the moment of inertia of the soliton. 

Using the expansions (|2l|) , (^2|) , in eqs. (|l9|), ( f20|) up to terms quadratic in R^R we arrive at a functional integral 
over the time dependent orientation matrices R(x4) with the action quadratic in angular velocities Q. In large _/V c 
limit, the latter corresponds to the hamiltonian of the quantum spherical rotator: 



H rot = J a 2 /(29) , 



(23) 



It means that despite of the fact that in general the path integral over R runs over all possible trajectories, in large 
N c limit the main contribution comes from trajectories close to those of the quantum rotator with hamiltonian (|23j). 
In eq. (pa) J a is the spin operator of the nucleon and 



9 = 9 s 



e 1 



is the total moment of inertia, including also the valence quark contribution 
to eq. (E3) the quantization rule (for the Euclidean time) is given by 



Sp(i?ti?r a ) 



Ja 

6 ' 



We also make use of the identity 



R?R=^Sp(R?RT a )T a 



(24) 

see eq. (|o|) below. According 

(25) 

(26) 



in order to separate the part, acting on the rotational functions of the nucleon iPn(R), from the Pauli matrices r a 
which act in the isospin space of the quark states. 

It should be also noted that in expansions (0) and (H we have neglected the terms containing derivatives of 
angular velocity <9 T (i?^i?). Such a treatment is consistent in a sense that for the hamiltonian (p^) we have 



d T (R^R) = [H rot ,R^R] = 



(27) 
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for any rotation state of H rot . 

After quantization ( p5|) we have to deal with not-commuting collective operators. Apparently, since the path integral 
corresponds to the following operator construction: 



R{T 2 )=R 2 



J VR Fi(i2(ti)) ■ • • F n (R(t n )) e~ s ^ = (R 2 ,T 2 | T{Fi(JZ(ii)) • . . F n (R(t n ))} \ ife, Ti) , 



(28) 



R(T 1 )=R 1 



where T stands for time-ordering of operators F n (R(t n )) in the Heisenberg representation, the order of the collective 
operators, which appear after the integration over R, is strictly fixed by the time ordering. 

Back to the sea contribution ([l9]), expanding both the action and the integrand in powers of R^R and performing 
the functional integration over time dependent R{xi) in large N c limit, we now have 



(N',p'\¥o^\N,p) sea = N c / d 3 xe^-^ / dR^* N ,(R)Sp\R^OR(x ) 0\ 



D(U) + WR 



,0) 



Sp(fitflr a ) j. j a /e 



(29) 



where ^ means again time-ordering of the collective operators. The rotational wave functions ipN^R) comes from 
the product of N c matrices R(T/2), which appear in eq. ( |l3|) after transformation (|l^), contracted by the matrices 
r^r from the definition of the baryon current (10|). For given spin J, J3 and isospin T, T3 these wave functions can be 
expressed through the Wigner D functions 



i>TT 3 j J3 (R) = (-i) T+Ts V2rTTD^j(R) . 



(30) 



In eq. (g9|) the integral over x has appeared due to the integration over the translational modes of the soliton and the 
matrix element in the integrand is given by 



Sp 



R^OR(0) (x,0\ 



-1 



D{U) + R*R 



x,0) 



Sp 



s P (-RtflT a ) j j a /e 



R ] OR (x,0| 



1 



D(U) 



\x,0) 



2 ^Sp J d A y[0{y A )J a R^OR + e{~ yi )tfORJ a } (f, 0| 7 ^|y)r Q (y| 7 ^ T |.r. 



>D(U)> 



<D(U)< 



As a next step we use the spectral representation for the quark propagator 

e„>0 



D(U) 



6{x' A - Xi) e- £ "«- x ^ $„(£') & n (x) - 6(x 4 - x' A ) £ e -^*-^) $„(£' ) $t (£) 



e n <0 



written in terms of eigenvalues e n and eigenfunctions <!>„ of the Dirac hamiltonian (0) 

h(U)$ n = £„$„ . 



(31) 



(32) 



(33) 



in order to evaluate the matrix elements and the integral over y 4 in (|3l]) . Finally we get for the sea contribution ( |l3| ) 
the following result: 

(N',j/\^0^\N,p) sea = N C /dV^ / dR^* N ,(R)\ V <(f)i? t Oi?$„(f) + -^-^^ 

^e„,<0 m,n 

x{-0(e n )0(-e m )J o [<(x)i?tOi?$ m (f)] +9K)«( £m ) R*OR$ m (x)] J a } x (to|t» j^C-ft) (34) 

Now the valence level contribution to the form factor (|l4|) is in order. Note that in the limit of large Euclidean 
time T — ► +00 the quark propagator (|32| ) is dominated by the contribution of the lowest level with e„ > 0, i.e. by the 
valence quark level with the wave function <S> va i(x*). We perform analogous calculations as in the case of Dirac sea, i.e. 
using the rotating ansatz ( |l6| ) in ( |l4| ) and expanding up to the linear order in the angular velocity with time-ordering 
of collective operators we arrive at the following result for the valence quark contribution: 



G 



(N'J\^CM\ N, p)™ 1 =N C J d 3 xc 1 ^ J dR i>* N > (R) | *J„i (*) ^ f Oi? (2) 

+^ E — — {^[^[^(^Oi^aK^M^ 

+0(-£„)[j a [<^(£)i?+OE $„(*)] (n|r>af> + [^ n (x) ^OR^ val (x)] (val\ T a \n) J a ]}^ N (R) (35) 
Finally we add the valence contribution (|3q ) to the sea part (|34|): 

(N\p'\^0^\N,p) = N c J d^e^-^ J dR^* N ,(R)\ ^ n {x) R^OR<P n (x) 

E {ja[$.J l (^)i? t Oi?$ m (f)](TO|T a |n) + [$] n (f)i?t 0i? $ n (f)]( n | ra | m ) jlL^fl) (36) 



It is easy to see that the result has the same structure as the sea contribution but with the valence quark included 
into the occupied states. 



III. REGULARIZATION 



The NJL model is not renormalizable. Even in the leading order of the saddle point approximation one faces 
ultraviolet divergences. The integrand in the r.h.s. of eq. ( |l3| ) contains fermionic determinant T)et[D(U)] and a 
matrix element: 



SpO(0|^-|0) = T ^Trlog[ J D(C7)-eO] 



e=0 



1 D{Uy 1 5e(0) 

both divergent. Apparently, in order to make it finite it is sufficient to regularize 

TrlogL> £ = Trlog[D([/) - eO] = Trlog[L»(C7) + i? f R - eR^OR] 



(37) 



(38) 



In fact, we are interested in a particular operator O, for which only the real part of TrlogZ? e is divergent. The 
corresponding regularized sea contribution to the matrix element (O) has the form 



(N',j/\^0^\N,p) s r e e a =N c / d 3 ^'^ / dRiP* N ,(R) 



Se(x, 0) 



Re(Trlog D e ) 



i/> N (R) , (39) 



Sp(fit_Rr Q ) J J a /e 



where we used a proper-time regularization 

Re(Tr log D e ) reg = i[Tr log(DjZ). 



e }\reg 



1 [ ^ Tr 
2 J u 



c -uDlD t _ e -uDlD 



(40) 



Hence 



5e(x) 



Re(TrlogD e ) re3 = -- Sp R\x A )OR( Xi ){x\e 



,-DtD/A 2 n -l 



D-'x) 



- Sp i?t(x 4 )Ot jR ( X4 )( x |( jD t ) -i e -c^/A 2 | 



Here D stands for 



D = D e=Q = d T + h(U) + R^R 



(41) 



(42) 
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Note that we assume e(x) to be real but the matrix O is not necessarily hermitean. 
Using the expansion 



D~ 



1 



d T + h d T + h d T + h 
up to the linear order in the angular velocity we get 

1 

A- 

1 



(43) 







+ i_ j dae^^-" 2 )/^ RlRe( 1 - a ^- h2 V A2 (d T + h) + 



(44) 



After quantization (|25J) of the time-ordered collective operators and using the spectral representation for the quark 
propagator ( |32|) we get for the regularized Dirac sea contribution (^9|)in the form 

(N',0\&O*\N,&™ = N C J d 3 xe^-^ 3 j dRr N '(R){^2^i(£n,v) rfOR$ n (x) 
+ Jq S [^2 +) (£m: R^ORZjx)} + R^£ ra ,£„, i?t OR $„ (f )] J a 

The coefficient r/ = 1 for hermitean matrix O and 77 = — 1 for anti- hermitean one: 

O f =T]0. 

The regularization functions are given by 

l+„ ?dc eh(^+4)/A 2 ] 
— , 



n\T a \m) Ujv(-R). (45) 



(46) 



(47) 



^(Ermenil) = 7£ (±) (e m ,e n ) - riTZ {± \e n ,e r , 



(48) 



OO OO 



1 



1 



1 

A 2 " 



00 —00 
1 

dae 





±i(u> — u') + (iu> + e n )(iuj' + e m ) 



A 2 



(49) 



IV. ELECTROMAGNETIC FORM FACTORS 



In this section we apply the formalism, developed above, to the case of the electromagnetic form factors. We pay 
a special attention to the ultraviolet regularization to preserve the electromagnetic gauge invariance. Following this 
regularization prescription we evaluate the nuclcon matrix elements of the electromagnetic current in the NJL model 
and relate them to the Sachs form factors. 

In the previous section we have presented a scheme for evaluating the regularized determinant with a small external 
source eO (^7|). In the case of the electromagnetic form factors this external source is an electromagnetic field A^: 



D(U)-eO -> D(U,QA IM )= l4 [ 7li (d li -iQA IJ ,)+MU^} 



(50) 



Here Q is the quark charge matrix 
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2/3 \ _ 1 
-1/3 I ~ 6 + 2 



(51) 



Our proper-time regularization in the Euclidean space ( |40| ) includes operator D. Therefore, in order to preserve 
the gauge invariance in ( |40| ) the Euclidean electromagnetic field should be real. 

The complex conjugate Euclidean Dirac matrices 7* are connected to 7 M by some unitary transformation V 



It is easy to show that 



where 



j^V^V- 1 ( M = 1,2, 3, 4, 5) 



(Vt 2 )D(U, QA^iVr 2 )' 1 = [D(U, Q'A^)} 



Q' 



T 

Y 



(52) 
(53) 

(54) 



and the asterisk stands for the complex (not hermitean!) conjugation. 

Comparing Q j5l| ) and Q' (^3) we see that their isovector part remains the same whereas the isoscalar component 
changes the sign. The identity (|53|) shows that the isoscalar electromagnetic form factors originate from the imaginary 
part of Tr log D(U, QA^) whereas the isovector form factors are generated by the real part of Tr log £>([/, QA^). 

Since the imaginary part of Tr log D(U, QA^) is ultraviolet finite, for the isoscalar form factors we use directly the 
non-regularized expression (p6|). In this case the integral over R becomes trivial and reduces to the orthogonality 
condition for the rotational wave functions or to the standard matrix elements of the spin operator J a 

(n,4,?\&iw^\T 3 ,J 3 ,$=N c 5 T , T3 J d 3 xe^-^|^ j3 &n(x)Ha^n(x) 



40 ^y^h Yl 



-n & m 



[$t(x)i 7 47 M * m (^)] {m\T a \n) + [$J ra (f)z 7 47 At $„(f)] (n\ T a \m) 



(55) 



In the case of the isovector current, however, we have to use the regularized expression ( pig ) for the sea contribution, 
which leads to: 

-40" Yl [^2 +) (e™> £n, v)Ja D 3b {R) + n { 2 - ] (e m , e n ,r))D 3b (R) J a ] [& m (£) i lAlll r b *„(£)] (n\T a \m)^ N (R) , (56) 



where 



Dab(R) = ~ S P (i?Vi?T b ) 



(57) 



Note that i^/Ci^J h is hermitean for fi — 1,2,3 and anti-hermitean for /i = 4 which means that the coefficient ij 
(defined in (E6h) has the following values: 



1 if /x = 1,2, 3 
-1 if/x = 4. 



(58) 



Now the rotational matrix elements in (56) can be easily computed using 



(T 3 -J' 3 \J a \T 3 J 3 ) = -5 T , T3 {T a )r 3 j 3 , 



(59) 



(T^4\D ab (R)\T 3 J 3 ) = 



(60) 
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From eq.(56) we find: 



129 



gXll^A ^ +) (e m ,£„,7?) '(£,„, +ieab c (r c )^,7 3 ^ 2 +) (e m ,£„,ry) - ft 2 ^e™, £„,/?) | 



x [& m {x) i~ia^ h $n(x)} (n\T a \m) 
Under the transformation Vt 2 (153) we have 



(Vr^iVr 2 )- 1 =7t = iL (Vr 2 )T a (Vr 2 )-' = -(r a )* = -(r Q f , 



(61) 



(62) 



(Vr 2 )h(Vr 2 



h* = h 1 



(63) 



where the superscript T stands for the transposition in both the matrix indices and the coordinate space. The last 
equation shows that the eigenfunctions of h can be chosen in such a way that 



(Vr 2 )<S> m {x) = <f>* m (x) 

and hence, the following relations hold: 

^i{x)-f 4 -f^ m (x)langlem\T a \n) = i] $ I t „(f)7 4 7 p $„(f)(n|r Q |m) , 

<S>i(x)j4JnT a <i> m (x){m\T a \n) = -i] ^ m {x)j4l^T a ^ n (x)(n\T a \m) 



(64) 

(65) 
(66) 



with values for rj are given in (|5^ 

Using these symmetry properties of matrix elements we can simplify expression (|55|) 



*|T 3l J 3 ,p) = N C 6 T , T3 J A^xe^-^hj, . h LJL ]T $t (2) z 74 7m 



(67) 



Starting from the general formula ( pq ) analogous calculations of the non-regularized valence quark contribution to 
the matrix element of the isovector part of the electromagnetic current can be performed: 

-i E I X _, [e(en)-vO(-e n )][(l-v)SabS A j 3 + {l + v)ieabc(T c )j^ 3 } 

x [$t (£) ^ 47aiT " $ oaJ (f)] ( wo J|r» I . (68) 



Up to now we worked with Euclidean Dirac matrices. Since the same combination i 7 47 M enters both left and right 
hand sides of our formulas for form factors the transition to the Minkowskian Dirac matrices is trivial — it is sufficient 
to replace the Euclidean i 747 ^ by Minkowskian ij ^^ everywhere in eqs. (|6l"|), ( |67[ ) and (|6^). Henceforward, we work 
in the Minkowski space-time. 

The nucleon electromagnetic Sachs form factors are related to the matrix element of the electromagnetic current 



f = * 7 M 0* = *t 7 <V ( - + -t 3 ) * 

6 2 



(69) 
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in the standard way: 



(4,P / \j (0)\J3,p) = G E (q 2 )5j,, h , 



(70) 



(J>'|j fe (0)|J 3 , P ) = 



2M 



N 



l (r l )ju 3 q m G M (q 2 ).. 



where 



q = p -p. 

The isoscalar and isovector parts of the form factors are defined by 

E(M) — 2 B(JVf) + J 3<J_E(M) • 



(71) 



(72) 



(73) 



Using our result (67) for the nucleon matrix elements of the isoscalar vector current and definition ( |70j ) we find for 
the isoscalar electric form factor a simple result: 



gFV) 



fd 3 x^ s { £ <(f)$„(x)- 



e n <e val e^'<0 

Here we have to subtract the vacuum contribution which provides the correct normalization: 

cr°(o) = i 

For the electric isovector form factor frorn(|6l"|) and ( |6q ) we obtain: 



(74) 



(75) 



E 



Nc_ 
68 

1 



where the regulator has the form 



^val &n 



' {x)r a $ va i{x)\ (val\T a \n) 



(76) 



1 



du I 1 e 



£" m 6 



(77) 



1/A 2 



It can be easily seen that the calculation of the moment of inertia is very similar to the case of the electric isovector 
form factor. Indeed, the Dirac sea contribution to the moment of inertia (Oil) can be written in the form: 



Q sea 6„. h / dr = 



N c d 2 
4 de a de b 



ReTrlog(Dt£) e ) 



reg j 



where now 



D e = D(U) - ie c r c . 



(78) 



(79) 



Note that the first e derivative in (|78|) can be interpreted as the derivative of the effective action with respect to the 
source to obtain the electric isovector form factor (at q 2 = 0) whereas the second e derivative is analogous to the 
expansion in angular velocity. A detailed calculation including the valence quark contribution gives: 



Q5 ab = ^\Y,{Mr a \n){n\r h \m)'R^{e mi e n )+ £ — ^ 

^ I , , , ^val 



[n\T a \val)(val\T b \n) 



(80) 
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This regularized expression has been already derived in ref. |7J and numerically evaluated in refs. |8Ul 
One can see from ( |76|) and (|8p| ) that the isovector electric form factor has a proper normalization : 



cr l (o) = i 



(81) 



as it should be. 

Now we proceed to compute the magnetic form factors. From matrix element ( |67| ) and the definition of the magnetic 
Sachs form factor (|7l]) we find that the first non-vanishing contributions come from rotational corrections: 



T =o (n 2^ _ N C M N kaj iq_ 
69 \q 2 



d 3 xe* 9 " J2 



$tn(xh l k ®n{x) (n\T a \m) 



(82) 



e m <e„„, 



- m c n 



The calculations of the magnetic isovector form factor are more involved and the final expression includes leading 
order terms as well as rotational corrections (next to leading order): 



T=1 2 N C M N kbj %q 
g m W ) = 3 ' — 



29 



£abc 



x) nr \m) 



£ £li^I[$t (f)7 o 7 V«$ w (f)]Mb 



(83) 



Here, two different regularization functions appear: 



^Ml( £ n) = -K-i{e n , 1) = £ n / — 



dw e -(^+e 2 J/A 2 



2vr w 2 + e 2 J y/u 

l/A 2 



(84) 



_ 1 f 1 d{3 (1- P)e m - (3e n n _ l(1 _ /3)£ 2 n+f3s 2 J/A 2 



4tt y ^(1 - /?) (1 - /JJe*. + /fe 2 



(85) 



Eqs. (|74|), (f76|), ( jS2| ) and ( |83] ) are our final expressions for the electromagnetic isoscalar and isovector form factors 
in the NJL model. According to ( |Z3| ) the proton and neutron form factors are expressed in terms of the isoscalar and 
isovector form factors as follows 



r v _ 1 \r T=0 _i_ r T=l l 

E{M) — Til- + E(M)\ J 



(86) 



(87) 



V. NUMERICAL RESULTS 



In the numerical calculations we use the method of Ripka and Kahana |16[ for solving the eigenvalue problem in 
finite quasi-discrete basis. We consider a spherical box of large radius D and the basis is made discrete by imposing 
a boundary condition at r = D. Also, it is made finite by restricting momenta of the basis states to be smaller than 
the numerical cut-off K max . Both quantities have no physical meaning and the results should not depend on them. 
The typical values used are D ~ 20/M and K max ~ 7M. 

In addition, all checks concerning the numerical stability of the solution with respect to varying box size and choice 
of the numerical cut-off have been done and the actual calculation is completely under control. 

The parameters of the model are fixed in meson sector in the well known way || to have f n — 93 MeV and 
m-n = 139.6 MeV. This leaves the constituent quark mass M as the only free parameter. 

The proton and neutron electric and magnetic form factors are displayed in Figs.||j^. The theoretical curves resulting 
from the model are given for four different values of the constituent quark mass, namely 370, 400, 420 and 450 MeV. The 
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magnetic form factors are normalized to the experimental values of the corresponding magnetic moments at q 2 = 0. 
With one exception of the neutron electric form factor (Fig.^), all other form factors agree with the experimental data 
quite well. The best fit is for the constituent quark mass around 420 MeV. 

As can be seen the only form factor which deviates from the experimental data is the neutron electric form factor 
and this requires some explanation. Obviously, this form factor is the most sensitive for numerical errors. According 
to the formula (|87]) the form factor has been calculated as a difference of the electric isoscalar and electric isovector 
form factors. Both form factors were of order of one and calculated by the code with high enough accuracy. However, 
the resulting neutron form factor has experimental values of order 0.04, i.e. about 4% of the value of its components. 
This means that even small theoretical uncertainties for one of the components can be enhanced by a factor 50. It 
means that the numerical accuracy together with the used large 1 /N c approximation behind the model are strongly 
magnified and could in principle result in a deviation from the experimental data for momentum transfers above 100 
MeV. It should be also stressed that the experimental data (see ref. |ft8| ), available for the neutron electric form 
factor, are strongly model dependent and a different N - N potential used in the analysis of the data can lead to an 
enhancement of the experimental numbers by more than 50 %. There are also very recent experiment Jl9[] , which 
deviate from the previous data (is|l , indicating larger values closer to our numbers. 

As the next step, we compute other electromagnetic observables: the mean squared radii, the magnetic moments 
and the nucleon-A splitting. In particular, the charge radii and the magnetic moments can be obtained from the 
form factors: 



(r 2 ) T=0 ,i 



6 dG 



'T=0,1 



G 



T=0,1 
E 



dq 2 



(88) 



^ = G T =°-\ q 2 ) 2 • (89) 

For the quark masses 370,420 and 450 MeV and pion mass m, = 140 MeV the calculated values are presented in 
Table Q. The values for the axial vector coupling constant 17,4, taken from ref. pT| , are also presented for completeness. 

The results of Table | (m n = 139.6 MeV) again indicate the value ~ 420 MeV for the constituent quark mass, in 
agreement with the conclusion drawn from the form factor curves. The same value has been suggested earlier P,[l0|], 
where a smaller number of observables has been considered. With the exception of the neutron electric squared radii, 
to which remarks similar to the case of the neutron electric form factor are valid, the contribution of the valence 
quarks is dominant. However, the contribution of the Dirac sea is non-negligible and it varies within the range 15 - 
40%. As can be seen, the numerical results for the nucleon N-A mass splitting (Ma — Mm), the mean squared proton, 
isoscalar and isovector electric radii and the axial coupling constant (gA) as well as the g-dependence of the proton 
electric and magnetic, and neutron magnetic form factor differ from the experimental data by no more then about 
±5%. Finally, for the magnetic moments we have got results 20-25% below their experimental values. However, 
despite of this underestimation of both magnetic moments, we have found a very good result for the ratio fi p / fj, n 
which is far better than in other models. The experimental ratio is almost exactly reproduced for the constituent 
quark mass 420 MeV. 

In the present calculations, in contrast to the isoscalar magnetic moment the isovector one includes non-zero 
contributions in both leading (N°) and next to leading order (l/N c rotational corrections). The enhancement due 
to these corrections improves considerably (see also ref. [^1)) the agreement with the experiment and resolves to a 
great extent the problem of strong underestimation p| of both the isovector magnetic moment and the axial vector 
coupling constant in leading order. It is also important to note that since the 1/N C rotational corrections have the same 
spin-flavor structure like the leading term they do not violate the consistency condition of Dashen and Manohar ]2Cj ] 
derived in the large- N c of QCD. The latter means that other 1/N C corrections (e.g. meson loops) should contribute 
to both the leading and next to leading order. 

In Table [n] we give the theoretical values of the same quantities but with the physical pion mass set to zero. The 
chiral limit (m w — ► 0) mostly influences the isovector charge radius. In fact, as it should be expected | pl| , p2[ , the 
isovector charge radius diverges in chiral limit and our calculations with zero pion mass confirms it. It can be seen 
in Fig.|6|, where the isovector charge radius is plotted vs. the box size D. The Dirac sea contribution to the radius 
grows linearly with D and diverges as D — > 00. Because of this that quantity (and the relative quantities) is not 
included in Table [n| The other observable strongly influenced by the chiral limit is the neutron electric form factor. 
For the m v — > the discrepancy from the experiment is by almost a factor two larger than in the case ^ 0. The 
other observables differ in the chiral limit by about 30%. The comparison of the values in the two tables indicate 
that taking the physical pion mass gives us a best fit with a much better agreement with the experimental data. In 
addition, in chiral limit we observe much larger contribution from the sea effects, about 50% of the total value. 
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In section 3 we have found that the magnetic isovector form factor is the only one which includes non-zero con- 
tributions in both leading (AT?) and next to leading order (1/N C ) rotational corrections. However, the numerical 
calculations show that the (1/N C ) rotational corrections do not affect the g 2 -dependence (slope) of the form factor 
but rather the value at the origin, G|^ 1 (0) = /i T=1 which is the isovector magnetic moment. It is not surprising 
since (as can be seen from eq. (^3|) in both leading and next to leading order terms the shape of the wave functions 
<!>„ determines the g 2 -dependence of the form factor whereas the value at q 2 — depends on the particular matrix 
elements included. In leading order the isovector magnetic moment is strongly underestimated Jl^| . As has been 
shown in ref. jll] the enhancement for this quantity due to the 1/N C corrections is of order (N c + 2)/N c and it 
improves considerably the agreement with experiment. However, as can be seen from Table [| the magnetic moments 
are still below the experimental value by 25 %. 

The isoscalar and isovector electric mean squared radii are shown in Figs.^j8| as functions of the constituent quark 
mass. The same plot for the proton and neutron electric charge radii is given in FigJ|. In these plots the valence 
and sea contributions are explicitly given (dashed and dash-dotted lines, respectively). As could be expected, for 
the isoscalar electric charge radius the valence part is dominant (about 85%), due to the fact that there is no 1/N C 
rotational contributions. This is not true for the isovector electric charge radius, where the sea part contributes to 
about 45% of the total value (Fig j^) . Also, this effect can be seen at the proton and neutron charge radii (Fig||) which 
are linear combinations of the isoscalar and isovector ones. For proton the sea contribution is about 30% whereas for 
the neutron charge radius the negative sea part is dominating and the valence contribution is negligible. 

In addition, in Fig.[l^ we plot the proton and neutron charge distribution for the constituent quark mass M = 420 
MeV. For the proton we have a positive definite charge distribution completely dominated by the valence contribution 
whereas in the case of the neutron the Dirac sea is dominant. In accordance with the well accepted phenomenological 
picture one realizes a positive core coming from the valence quarks and a long negative tail due to the polarization of 
the Dirac sea. Using the gradient expansion the latter can be expressed in terms of the dynamical pion field - pion 
cloud. 

For completeness, in Fig.|ll| we present also the magnetic moment density distribution for proton and neutron for 
the constituent quark mass M = 420 MeV. The sea contribution becomes non-negligible for distances greater than 
0.5 fm. Due to the relatively large tail contribution to the magnetic moments the sea contribution to these quantities 
is about 30%. 



VI. SUMMARY AND CONCLUSIONS 



Our numerical results support the view that the chiral quark soliton model of the SU(2) Nambu-Jona-Lasinio type 
offers relatively simple but quite successful description of some low-energy QCD phenomena and, in particular, of the 
nucleon electromagnetic properties. In calculations, using f v = 93 MeV, = 139.6 MeV and a constituent quark 
mass of M = 420 MeV, and no other parameters, we have obtained an overall good description for the electromagnetic 
form factors, the mean squared radii, the magnetic moments, and the nucleon-A splitting. The isoscalar and isovector 
electric radii are reproduced within 15%. The magnetic moments are underestimated by about 25%. However, their 
ratio Hp //J, n is almost perfectly reproduced. To this end, the account of the 1/N C rotational corrections is decisive. 
The g-dependence of G E (q 2 ), G p M (q 2 ) and G^ t (q 2 ) is very well reproduced up to momentum transfer of 1 GeV. The 
only exception is the neutron form factor G E {q 2 ) which is by a factor of two too large for q 2 > 100 MeV 2 . One should 
note, however, that G E (q ) is more than an order of magnitude smaller than G E (q 2 ) and as such it is extremely 
sensitive to both the model approximations and the numerical methods used. Here, it turns out that the agreement is 
noticeably worse if the chiral limit (m n — > 0) is used. This can be easily understood since the pion mass determines 
the asymptotic behavior of the pion field. Altogether we conclude that the chiral quark soliton model based on a 
bosonized SU(2) NJL type lagrangean is quite appropriate for the evaluation of nucleonic electromagnetic properties. 
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FIG. 1. Diagrams corresponding to the expansion in Q of the current matrix element: a) the valence contribution and b) 
the Dirac sea contribution. 

FIG. 2. The proton electric form factor for the momentum transfers below 1 GeV. Experimental data from ref. Jl7[ |. 



FIG. 3. The neutron electric form factor for the momentum transfers below 1 GeV. experimental data from refs. |18j,|l 

FIG. 4. The proton magnetic form factor normalized to the experimental value of the proton magnetic moment at q 2 = 
for the momentum transfers below 1 GeV. The normalization factor can be extracted from Table 1. Experimental data from 
ref. @. 

FIG. 5. The neutron magnetic form factor normalized to the experimental value of the neutron magnetic moment at q 2 = 
for the momentum transfers below 1 GeV. The normalization factor can be extracted from Table 1. Experimental data from 
ref. @. 

FIG. 6. The isovector charge radius as a function of the box size D for m w = and = 139.6 MeV. 

FIG. 7. The isoscalar electric charge radius as a function of the constituent quark mass M. The valence and sea parts are 
marked by the dashed and dashed-dotted lines. 

FIG. 8. The isovector electric charge radius as a function of the constituent quark mass M. The valence and sea parts are 
marked by the dashed and dashed-dotted lines. 

FIG. 9. The electric charge radii of proton and neutron as functions of the constituent quark mass M. The valence and sea 
parts are marked by the dashed lines. 

FIG. 10. The charge density distribution of the proton (lower) and neutron (upper) for the constituent quark mass M = 420 
MeV. 



FIG. 11. The magnetic moment density of proton and neutron for the constituent quark mass M — 420 MeV. 
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TABLE I. Nucleon observables calculated with the physical pion mass. 



Quantity 


Co 

370 MeV 


nstituent Quark 1VI 

420 MeV 


ass 

450 MeV 


Exper. 


total 


sea 


total 


sea 


total 


sea 


< r 2 >t=o [fni 2 ] 


0.63 


0.05 


0.52 


0.07 


0.48 


0.09 


U.DZ 


< r >t=i [im J 


1.07 


0.33 


0.89 


0.41 


0.84 


0.45 


U.OD 


^ 2 ^ r-f 21 


0.85 


0.19 


0.70 


0.24 


0.66 


0.27 


A T A 

0.74 




-0.22 


-0.14 


-0.18 


-0.17 


-0.18 


-0.18 


0.12 


r l 

^ T =o [n.m.J 


0.68 


0.09 


0.62 


0.03 


0.59 


0.05 


0.88 


r l 

[i T=1 [n.m.J 


3.56 


0.77 


3.44 


0.97 


3.16 


0.80 


4.71 


^ip [n.m.] 


2.12 


0.43 


2.03 


0.50 


1.86 


0.43 


Z. t\) 


Mn [n.m.] 


-1.44 


-0.34 


-1.41 


-0.47 


-1.29 


-0.38 


— i.yi 


i . . / , . i 


1.47 




1.44 




1.44 




1 a a 
1.46 




1.08 


0.32 


0.66 


0.28 


0.56 


0.25 


U. /4 


< r 2 [fm 2 ] 


1.17 


0.51 


0.65 


-0.31 


-0.52 


-0.24 


-0.77 


M A - M N [MeV] 


213 




280 




314 




294 




1.26 


0.08 


1.21 


0.11 


1.13 


0.06 


1.26 


TABLE II. Nucleon observables calculated with the zero pion mass. 


Quantity 


Co 

370 MeV 


nstituent Quark 1VI 

420 MeV 


ass 

450 MeV 


Exper. 


total 


sea 


total 


sea 


total 


sea 


< r 2 > T=0 [fm 2 ] 


0.88 


0.20 


0.66 


0.26 


0.61 


0.23 


0.62 


^T=o [n.m.] 


0.66 


0.07 


0.59 


0.09 


0.57 


0.09 


0.88 


^T=i [n.m.] 


4.61 


1.59 


4.38 


1.89 


3.91 


1.48 


4.71 


Hp [n.m.] 


2.63 


0.83 


2.49 


0.99 


2.24 


0.79 


2.79 


M« [n-m.] 


-1.97 


-0.76 


-1.89 


-0.90 


-1.67 


-0.70 


-1.91 


\Vp/Vn\ 


1.34 




1.34 




1.34 




1.46 


M A - M N [MeV] 


221 




261 




301 




294 
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G^(Q 2 ) [n.m 




sovector magnetic moment (n.m.) 




G) 

o 
o 




R [fm] 



